Superconducting fluctuations in granular metals with a large coupling between the 

grains 
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We study the fluctuation conductivity of superconducting granular metals at low temperatures 
and strong magnetic field destroying the Cooper pairs. Explicit calculations are performed for larger 
values of the coupling between the grains than those considered in previous works. We show that in 
a broad region of the coupling constants the superconducting fiuctuations still significantly reduce 
the conductivity leading to a negative magnetoresistance. 
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I. INTRODUCTION 



During the last decade study of electric properties of 
non-homogeneous metals attracted a lot of attention. In 
particular, granular metals have been investigated re- 
cently in a number of experimental worksldij. Theory 
of superconducting fluctuations in the graruilated super- 
conductors was suggested recently in Ref J3a. 

In these works a three-dimensional (3D) array of grains 
placed in a strong magnetic field H > He, where He is the 
field destroying the superconducting gap in the grains, 
and at low temperatures T <ti Te, where Te is the super- 
conducting transition temperature, was considered. It 
was demonstrated that the correction to the fluctuation 
conductivity is negative and this effect should exist even 
at zero temperature. The resistivity increases and only 
at extremely strong magnetic fields, H ^ He, reaches 
its classical value. Therefore the system exhibits a nega- 
tive magnetoresistance. This theory may explain existing 
experimental. 

It was important for the calculations presented in Ref.EI 
that the dimensionless conductance satisfied the condi- 
tion g ^ Ao/(5, here Aq is the BCS gap at zero magnetic 
field and S is the mean level spacing. 

The same effect of the negative magnetoresistance was 
obtained in a recent paperB for two-dimensional (2D) 
homogeneous superconducting samples at low temper- 
atures, T ^ Te, and strong magnetic field, H > He- 
The limit of the homogeneous metal is opposite to the 
one considered in the worksQ'cl because the dimensionless 
conductance g of a homogeneous sample is proportional 
to kpl, where kp and I are the Fermi momentum and 
the mean free path, respectively, and can be very large. 
The negative magnetoresistance can also be seen-uader 
certain circumstances in High-Tc superconductorsflij for 
high temperatures, T > Te, and low magnetic fields, 

H < He. nn 

In the present paper we generalize the results of Ref.l3'cl 
to larger values of the tunneling dimensionless conduc- 



tance g. In particular, the assumption that the dimen- 
sionless conductance, g, is restricted from above by Ao/(5 
is now dropped. This means that the structure of the 
granular metal becomes more similar to that of a bulk 
metal. The main question we are dealing with in this 
paper is if the superconducting fluctuations may cause 
a negative magnetoresistance in the granulated systems 
with the larger coupling between the grains. At the same 
time, our region of parameters is different from the limit 
of a homogeneous metal. We assume that 



1< .g < Et/S, 



(1.1) 



where Et = DjB? is the Thouless energy of a single 
grain, D is the diffusion coefficient and R is the radius 
of the grain. The last inequality in Eq. (1.1)^, although 
being more general than that used previouslyou, means 
that the granular structure is still important for our con- 
sideration. 

The granular material we consider now consists of a 
3D array of metallic grains with a typical diameter of 
the grains of 100 ± 20A. The electrons can tunnel from 
one grain to another. It is this tunneling that determines 
the properties of the entire system. Inside the grains 
there can be impurities and the shape of each grain is 
not perfect, so that the electrons are scattered randomly 
by the boundaries. Since the hopping amplitude is not 
very large, the macroscopic charge transfer is determined 
by the ratio of the hopping amplitude t to the mean level 
spacing b, or, in other words, by the dimensionless con- 
ductance 3 = X (I) ■ '^'^ limit t ^ S the discreteness 
of the energy spectrum in a single grain is not resolved 
and therefore the electron motion is diffusive through 
many grains. This limit corresponds to a macroscopi- 
cally weak disorder and results in a large dimensionless 
conductance g ^ 1. Below, we restrict our consideration 
by this limit. 

Let us discuss what happens with a granular metal at 
low temperatures. Below the critical temperature, Te, 
the electron-phonon interaction leads to the formation of 
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a superconducting gap in each grain and Cooper pairs 
appear. Applying a strong magnetic field one destroys 
the superconducting gap in each grain and comes to the 
picture of a normal metal with superconducting fluctua- 
tions. Our calculations are performed in this regime. 

We assume that the energy parameters are ordered as 
follows 



(1.2) 



The last inequality in Eq. (1^) means that the size of 
a single grain, R, is much smaller than the coherence 
length ^0- In this limit the superconducting fluctuations 
in a single grain are zero dimensional. We want to em- 
phasize that all energies are smaller than the Thouless 
energy, Et, and as a consequence the behavior of the 
system does not depend on grain boundaries or on indi- 
vidual scattering processes. Due to the large values of 
the conductance 5 3> 1 we may neglect weak localization 
and charging effects. Therefore all the effects considered 
below are entirely due to the superconducting fluctua- 
tions. 

The superconducting pairing inside the grains can be 
destroyed by both the orbital mechanism and the Zee- 
man splitting. The critical magnetic field destroying 
the superconductivity in a single grain in this case can 
be estimated as H°^R^ sa (f)Q, where (po = he/ e is a flux 
quantum, R is the radius of a single grain and ^ = a/^qI is 
the superconducting coherence length. The Zeeman crit- 
ical magnetic field can be written as g^sH^ — Aq, 
where /i^ is Bohr's magneton and g is the Lande factor. 
The ratio of this two fields can be written in the form 
H°'/H^ « RJR, where R^ = ^(poO"^- For R > Rc 
the orbital critical magnetic field is smaller than the Zee- 
man critical magnetic field < and the super- 
conductivity is suppressed by the orbital motion of elec- 
trons. Although the Zeeman mechanism can be easily 
included in the present consideration, we consider now 
only the orbital mechanism of the destruction of the su- 
perconductiwity. This limit is opposite to the one consid- 
ered in Ref.cl, where the Zeeman splitting was assumed 
to be the main mechanism of destruction of the Cooper 
pairs. A broader region of the conductance g used in 
the present paper makes the calculation somewhat more 
difficult than previously because one has to consider ad- 
ditional diagrams and calculate them using more compli- 
cated expressions for integrands. The remainder of the 
paper is organized as follows. In Sec. |l| we formulate the 
model. In Sec. [II we discuss the fluctuation conductivity 
of granular metals. In Sec. IV we discuss the Weak Lo- 
calization correction to conductivity of granular metals. 
Our results are summarized in the conclusion. 



II. THE MODEL 

We assume that the grains are packed in a 3D lattice 
surrounded by an isulator. The grains are coupled with 



each other and therefore the electrons can hop from one 
grain to another. Inside the grains as well on the surface 
there can be impurities and the electrons can interact 
with phonons. The Hamiltonian of the system can be 
written in the form 



H — Hq + Ht, 



(2.1) 



here Hq describes a single grain with electron-phonon in- 
teraction in the presence of a strong magnetic field and 
is given by 



Hq — ^ i?i,fca| j,ai,fc — |A| ^ a.i,k'^i,-k'^i,-k'ai,k' + Hi^p, 



(2.2) 



where i stands for the number of the grain, k = (k, "f 
),— A; = (— k, J,). The quantity A is an interaction con- 
stant and Himp describes the elastic interaction of the 



electrons with impurities. The interaction in Eq. (2.2) 
contains diagonal matrix elements only. This form of 
the interaction can be used provided the superconduct- 
ing gap Aq satisfies the last inequality in Eq. (1.2). The 
second term in Eq. (2T) describes tunneling of electrons 
from grain to grain and is given by 



-I- C.C, 



(2.3) 



where A is the vector potential, dy is a vector connect- 
ing the center of grain i with the center of a neighboring 
grain j (|dy| = 2R). The operator aj^ is the creation- 
operator of an electron in grain i in the state p and Uip is 
the annihilation operator of an electron in grain i in the 
state p. 



III. FLUCTUATION CONDUCTIVITY 

In this section we consider the conductivity of granular 
metals in detail. The dc conductivity a is ijeliated to the 
operator of the electromagnetic response aai3: 



a = Imi , 

w^o —iio 



(3.1) 



where Q^'iuS) is the analytical continuation of QiiiOy) 
into the upper complex half plane and is called the re- 
tarded operator of the electromagnetic response, lo is 
the frequency of the external electromagnetic field. In 
order to calculate Qiiijjy) we use Matsubara's diagram 
techniqueEj. After calculation of Qiiuiy) for imaginary 
frequencies we have to carry out the analytical con- 
tinuation of Q(iiOy) into the region of real frequencies: 
iu^y uj ^ iO^. All diagrams which contribute to the 
conductivity of the granular metal are shown in Fig. |^. 
The same clasSr^tdiagrams describe the conductivity of 
the bulk metal.ElO Scattering of the electrons inside the 
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grains by impurities is included in the Born approxima- 
tion, giving rise to a scattering mean free time t and 
resulting in a renormalization of the single electron nor- 
mal state Green's function to G°(ie„, p) = (ie„ — ^(p) + 
i/2Tsign(e„))~^, here e„ — {2n + 1)ttT is the fcrmion 
frequency and ^(p) = e(p) — Ef is the electron energy 
counted from the Fermi level. For / <C ic, where / is the 
mean free path and Lc is the cyclotron radius, we can 
treat the Green's function in the quasiclassical approxi- 
mation. In this approximation the magnetic field results 
in the appearance of an additional phase: 



X(ir!fc,q) = 



In 



£o{H) + \nk\ 



An 



(3.3) 



here fifc = 2knT is the boson frequency, vq is 
the density of states on the Fermi surface, /^(q) = 
8/37r((;J/Ao) ^-^^(1 — cos{qid)) describes the tunnel- 
ing of electrons from grain to grain; £o{H) = 
(2/5)(0/(/)o)^i?T, where (j) is the magnetic flux through 
the grain. Th e propagator of superconducting fluctua- 
tions, Eq. ( ^.3| ), is presented by the sum of all diagrams 
with two incoming and two outgoing lines in Fig. 



G{ien, r - r') = G°(ie„, r - r') exp 




Ads 



(3.2) 



Each wavy line in the diagrams represents the propa- 
gator of the superconducting fluctuations K{i^lk,'i)- 



/ 



+ 




\ 



FIG. 1. The Dyson equation in the ladder approximation 
for the propagator of the superconducting fluctuations. The 
black point represents the coupling constant and the shaded 
three-point vertex stands for the renormalized impurity ver- 
tex. 



J 




9 10 

FIG. 2. Diagrams for the leading order contribution to the fluctuation conductivity of granular metals. Wavy lines symbolize 
the propagator of the superconducting fluctuations, thin solid lines with arrows are the normal state Green's functions averaged 
over impurity positions and shaded semicircles are vertex corrections arising from impurities. Dashed lines with central crosses 
are additional impurity renormalizations and shaded blocks are impurity ladders. Diagram 1 is the Aslamazov-Larkin (AL) 
contribution, diagram 2 is the Maki-Thompson (MT), 5, 6, 7 and 8 are the density of states (DOS) diagrams. Diagrams 3,4 
and 9,10 arise when one averages the DOS and MT diagrams over impurities. 
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The impurity vertex entering these diagijatps is presented 
as a shaded half circle and has the formEla 



A. Correction to the Conductivity due to 
Suppression of DOS 



1 ei-eni^k - Sn)) 



(3.4) 



where 0(x) is the Heaviside- function and the third term in 
the denominator describes tunneling processes from grain 
to grain. Each external vertex is given by — 2etd sin(pd), 
where e is the charge of an electron, t is the tunneling 
amplitude and d is a vector connecting the centers of 
two neighboring grains. In the following subsections we 
consider the contributions to the fluctuation conductivity 
that arise from these diagrams. 



In this section we consider the correction to conduc- 
tivity due to suppression of the density of states (s. di- 
agrams 5 - 10). This (DOS) contribution arises from 
corrections to the density of states due to the super- 
conducting fluctuations. Even at strong magnetic fields 
{H > He) there are still some electrons that form fluc- 
tuational Cooper pairs. These electrons are bound and 
cannot simultaneously take part in one-electron charge 
transfer. This results in a reduction of the number of 
carriers for the one-electron charge transfer and the con- 
ductivity decreases. The analytical expression for the 
operator of the electromagnetic response taking into ac- 
count summation over the spin indices has the form: 



8 f 

Q(zw,) = -^r^ / — ^i^(ir!fc,q)T^C72(ie„,zrJfe-ie„,q)/(zr!fc,i£„+.), 



(3.5) 



where a factor of 2 originates from a similar diagram conductivity. The analytical expression for /(ifi^;, ie„+y) 
shown in Fig. which gives the same contribution to in Eq. (p^) is given by 



I{ink,ien+u) = j sin^ {q[d) ^^^c^^'^^^y T ^ j ^^-|^^^G(ie„+^, pi)G'^(z£„, P2)G(irifc - ie„, P2). (3.6) 



I{iQ.k,ien+v) 



here pi and p2 denote the momenta of electrons inside different signs. In all other cases the result equals to zero, 
the grains, q' is the quasimomentum and £n+u = £n + ^v For 17^, £„, uj^ -^I/t we obtain for I(iQki i^n+v) 
Inside the Green's functions we may put f2fe equal to 
zero because the characteristic frequency of the super- 
conducting fluctuation propagator, K{i^lk,qi)^ is of the 
order Jlfc ~ Aq which is much smaller than the Thou- 
less energy Et- The integral over p2 in Eq. ( ^.6[) is only 
nonzero when the poles of the Green's functions corre- 
sponding to one grain lie on different sides of the axis of Now usi ng Eq 
the real numbers. Therefore e„ and fife — £„ must have Eq- (p. 61) 



-enen+u)Q{ — £n(P'k — S-n))- 

(3.7) 



7f) we calculate the sum over £„ in 



D{iVLk,iuju,ai) = t'^C'^ {iVtk,i^k - «£„, q)/(if2fe, i£„+,y) 



-Eni^k - £«)) 



(|2£„-f!fc|-t-47rTaq)2' 



(3.8) 



Carrying out the summation over the Sn in Eq. (3^) we obtain for the function D{iQk^ i^u-, ^) the following result 



06 

D(jfifc, iLV^, q) = -^y— 6'(f7fe + uji,) 



47rT 



2 47rT 



(3.9) 
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where ip{x) is the logarithmic derivative of the Gamma- 
function and ttq is given by 



1 



AttT 



16 



(^(l-cos(g^d)) . (3.10) 



The second term in Eq. (3.1C) arises due to the renor- 
mahzation of the Cooperons when tunnehng processes 
from gra in to grai n ar e taken into account. We insert 
Eq. (3^) into Eq. (3^) and present the operator of the 
electromagnetic response in the following form: 



i=l Uk ^ ^' 

(3.11) 

Now the function Q{iuju) must be continued analytically 
into the upper complex half-plane of the frequency. The 
analytical continuation in Eq. ( 3.11 ) is carried out in Ap- 
pendix As a result, we obtain for the operator of the 
electromagnetic response Q{iuj^) after analytical contin- 
uation the following expression 



TT^dT^Aovo 



An ^ 



E 



(2^ 



coth 



2T ^ 

A? 



^y2(q) 



sinh^ -0- Bl 
sum 2T 



ry2(q) 2T 



, (3.12) 



where fir. 



Aq is an upper cut-off, ?7(q) — ri{q) + 2h and h 



Ha 



is the reduced magnetic field. In Eq. ( 3.12 ) 



we may put q = inside i^i^) because the main contribution to the integral over q comes from small momentum and 
in this case aq is a slowly varying function of q. The remaining integrals over Q in Eq. ( 3.12| ) can be easily calculated 
and the final result for the electromagnetic response is 



AnTj^J (2^)3 



In 



1 

2? 



^(0 + ev^'(c)-i 



(3.13) 



here we introduced the dimensionless parameter ^ = 
Aofi/AirT. For very low temperature T <C Aof], (^ » 1) 
we may use the asymptotic expansion for ip{^) and finally 
obtain the correction to the conductivity due to the sup- 
pression of the density of states: 



^DOS 



TT Aq 



1 



f7(q) 



(3.14) 



where 



(1;;^, and (To = is the clas- 



sical conductivity of a granular metal. One can see that 
the DOS diagram gives a negative contribution to the 
conductivity. The absolute value of is a decreasing 

function of the magnetic field H and reaches its maxi- 
mum value at the critical magnetic field H = He- lie 
absolute value of this maximum can be estimated asaQ 



^DOS 



(3.15) 



This maximum value is smaller than unity and this fact 
ensures our diagrammatic expansion. The conductivity 

^DOS is 

independent of the temperature and therefore 
it remains finite in the limit T ^ 0. This fact indicates 
that there are still virtual Cooper pairs even at zero tem- 
perature and strong magnetic field. The quantity a^^^ 
becomes comparable with ctq when g is of the order of 
unity. Such values of g mean that we would not be far 
from the metal-insulator transition. In this case we have 



to take into account all localization effects and Eq. (3.14) 
can be used only for rough estimates. 

We would like to note thai Eq. (3.14) does not differ 
from those written in Refs.Ellj. However, other contribu- 
tions to the conductivity considered in the next subsec- 
tions may change, so the extension of the calculations to 
the entire region specified by Eq. (1.1) is not as simple. 




FIG. 3. Higher order correction to the DOS 

As it was shown in Ref.BB we can neglect the higher 
order corrections to the DOS, an example is shown in 
Fig. H. At low temperatures T ^ Aory this diagram con- 
tains an additional small factor of (S/Aq) Infj. 

In the following sections we discuss the Aslamazov- 
Larkin (AL) and Maki-Thompson (MT) contributions 
to the fluctuation conductivity. It turns out that the 
AL contribution is proportional to whereas the MT 
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contribution can be divided into two parts. One is pro- B. Aslamazov-Larkin Contribution to Conductivity 

portional to at low temperatures and the other one 

remains finite when T 0. This contribution originates from the abijiLy of virtual 

Cooper pairs to carry an electrical currenill3. The dia- 
gram for the operator of the electromagnetic response is 
represented by the first diagram in Fig. ^ and its analyt- 
ical expression has the following form 



Q^^(ic^,) = ^^T^ [ ^^K{ink,<l)Kiin„-iu;,,c{)B^iink,iu;r..c{), (3.16) 

B{iil.k,ioJiy,q) describes the block of the Green's functions: 
B{ink,iu;,,q) = -i J j—^ sm{q^d) COs((g, ~ q^)d) ^^^^^^^ T ^ j -j^^ 

xG(i£„+^,Pi)G(irifc - ie„+^,pi)G(ie„,p2)G(ifJfe - i£„+^, p2)G(ie„+i,, ififc - ie„+i., q)G(ie„, ififc - iuj^^q), (3.17) 



here sin(gj'c?) denotes the current- vertex and cos((gi — 
ql)d) describes the tunneling-vertex. In this integral we 
may put flk = and uj,y = because in the vicinity of 
the critical magnetic field He the leading contribution to 
the response Q^^ arises from the fluctuation propagators 
rather than from the frequency dependence of the func- 
tion B, Eq. ( 3.17 ). Therefore we neglect the dependence 
of function B on frequencies flk and Lu,y. Integrating over 
the momenta pi,P2 and quasimomentum q', then sum- 
ming over the internal frequency e„ we obtain in the limit 
of low temperature for function -B(q) the result: 



s(q) 



. 32ged 



sin(g-d) cos{{q, - q[)d) 



d^q' 
(27r)3 



(3.18) 



As in Sec. [II A , in order to calculate the response Q 



AL 



we have to make the analytical continuation from the 
Matsubara frequency uii, to real values of uj. This c an be 
achieved by transforming the sum over Qk in Eq. (|3T|) 
into a contour integral. This procedure allows to carry 
out the analytical continuation of QiiiOy) into the upper 
complex half plane, ilOi, — > cj -|- iO"*". As a result for the 



sum over the Q,k in Eq. (3.16) we obtain 



QjTtT^ 1 

T^Kiink,ci)K{ink-iiu,,q) -^to-^ — . (3.19) 



Using Eq. (3.18) and Eq. ( |3.19 ) we finally obtain for the 
AL contribution to the fluctuation conductivity 



_ 64 S^T^ A 
^ 27^ ^ 

^ i—l 



sin^{qid) 



(3.20) 



which agrees with the result obtained in Refs.''u. One can 
see that a"^^ gives a positive contribution to the fluctua- 



tion conductivity and is proportional to T^, therefore it 
vanishes in the limit T — )■ 0. 



Let us estimate the right hand side in Eq. (3.20) for the 
case of low temperatures, T <Si Aq??, and near the crit- 
ical magnetic field, h < C qS /Aq. The main contribution 
to the integral in Eq. (3.20) comes from small momenta. 
Therefore we can make an expansion in q of siii{qid) and 
cos{qid). Retaining only the first nonvanishing terms and 
extending the range of integration from 27r/d to infinity 
we obtain 



-3/2_ 



.3/2, 



^1/2 \H-Hc 



3/2 



(3.21) 



The AL contribution grows when approaching the critical 
magnetic field He thus leading to a decrease of resistivity. 
In order to determine which contribution to the conduc- 
tivity will dominate we compare a"^^ with the maximum 
value of a^*^^ 



^AL 



rDOS 



9 



-3/2212 



Ay'53/2 



In" 



Ao 



Hr 



H-Hr 



3/2 



(3.22) 



For h <^ gS/Ao and sufficie ntly low temperatures, T <C 



Ao?y, one can see from Eq. (|3.22D that {a"^^ /a'^'^'=\ < 1 



This means that the AL contribution cannot change the 
monotonous increase of the resistivity of granular metals 
when decreasing the magnetic field. 



C. Maki-Thompson Contribution to Conductivity 

This contribution to the conductivity comes from co- 
herent electron scattering forming a Cooper pair on 
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impuritiesll2l,EJ. The MT contribution is represented by 
the second diagram in Fig. |^. The analytical expression 
for the operator of the electromagnetic response for the 
MT contribution to fluctuation conductivity is given by: 

i=l Ofc ^ ' 

(3.23) 



where K{iil.h,q) is the propagator of the superconduct- 
ing fluctuations and B(ink,iuj^,q) is a function describ- 
ing the contribution of the loop. This function has the 
form 



B{iVtk,i^u-: q) = / sin(gf,-d) sin((gfj - q[)d) 



(27r)3 {t:votY 



{2-Kf 



,_,,q). (3.24) 



In evaluating the sum over the Matsubara frequency e„ 
it is useful to break up the sum into two parts. In the 
first part e„ is in the domains ] — oo,—uj^[ and [0, oo[. 
This gives rise to the regular part of the MT diagram. 
The second (anomalous) part of the MT diagram arises 
from the summation over the e„ in the domain [—uJi,, 0[. 
Using this we can perform the sum over the £„ and after 



integration over the momenta we can write the function 
_B as a sum of an anomalous -B°" and a regular B^'^^ 
contribution to the MT diagram: 

2 

B{ink,ti^u, q) = -32^ cos(g,d)(B'^" -f B'^^'?), (3.25) 
where 



B" 



1 0{nk)0{io, - nk+i) 

47r Lu,, + AirTan 



^ ' 2 AnT 



^ ' 2 AttT 



(3.26) 



I 2Lo^ + n. \ fi ilk 



(3.27) 



Note that the anomalous and regular part have different 
signs and the anomalous part has an additional diffusion 
pole in comparison with the regular one. 



1. The Anomalous MT Contribution to Conductivity 



Q""(^c^.) 



E 



d^q cos(giC?)i^(t^,.,q) 



{2n) 



The function F^iiv^, q) is given bylll: 



(3.28) 



The term _B"" corresponds to the anomalous MT con- 
tribution, it appears in the case of a special p ole a rrange- 
ment in the integration over momenta in Eq. ( 3.24 ) , when 
the poles of the Green's functions corresponding to dif- 
ferent grains lie on different sides of the axis of the real 
numbers. In order to calculate the fluctuation conduc- 
tivity we have to insert into Eq. (|3.23 ). Then, one 
should perform the analytical continuation of Q°"(ia;,y) 
to real values of the external frequency uj. Finally, us- 
ing Eq. ( |3.lD we get the conductivity. To this end we 
represent Q"'" in the following form: 



fe, uji,, qj 



(3.29) 



haq)-t/'(^ + j^+aq) 



where f{ink,iuj^,q) o;,^i"/ci_, , , 

Bec ause of the presence of the Heaviside function 9 in 
Eq. ( |3.26| ), B"" is only nonzero in the domainjpi, 
The upper limit of the sum over fl^ in Eq. (3.29) depends 
on the external frequency a;^. Therefore it is not correct 
simply to make the s ubsti tution iuji, — > uj. Note that we 
have calculated Eq. ( 3.26| ) only for flk > but one can 
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easily obtain the corresponding expression for fl^ < re- 
placing f2fe by I rife I . Therefore, instead of summing over 
all values of fifc we extract the term ftk — and multiply 
the sum over fi^ > by a factor of 2. The prime on the 



sum in Eq. ( |3.29 ) indicates that the term with fife = 
should be multiplied by 1/2. The analytical continuation 
is achieved by transformation of the sum into Eueontour 
integral. The final result has the following formEJ: 



sinh 



2 _n_ 

2T 



+ a ) 
47rT ^ "qy 



(3.30) 



For ^ 1 we expand the numerator in a Taylor-series 
around = and confine ourselves to the first nonvan- 
ishing term. Then we obtain 



2iuj 



(47rr)2Ao 



Aq \ 



sinh' ^ ^ + 7)2 ■ 



(3.31) 



The last integral can be calculated in a sim ilar way as the 
second integral in Eq. ( 3.12| ) in Sec. Ill A and we obtain 
for the anomalous MT contribution to conductivity the 
following result: 



From Eq. (B.33) we see that (J^^ gives a positive contri- 



bution to the fluctuation conductivity and grows when 
approaching the critical magnetic field H He- The 
anomalous MT contribution is proportional to as tem- 
perature goes to zero. At zero temperature the anoma- 
lous MT contribution vanishes. 

It is interesting to compare cr*^"^ with the contribution 
to conductivity that arises from the suppression of the 
density of states a^'-'^ . The ratio of these two quantities 
is given by 



yDOS 
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90 



He 



H-H, 



1/2 



(3.34) 



9A2 E\4^raq,72(q) 



/ cos{qid) 



(3.32) 



Let us estimate the anomalous MT contribution to the 
conductivity. The main contribution to the integral in 
Eq. (3.32) comes from small momenta q, therefore for 
h <C gS / Aq we expand the numerator and denominator 
in powers of q. Confining ourselves to the first nonvan- 
ishing order in q, we obtain the following result: 



-3/2_ 



T2 



H-Hr 



1/2 



(3.33) 



One can see that at T ^ AqT? the anomalous MT con- 
tribution is small compared to <Jmax ■ Thus, we conclude 
that the anomalous MT contribution cannot change the 
monotonous increase of the resistivity of a granular su- 
perconductor when decreasing the magnetic field. 



2. The Regular MT Contribution to Conductivity 

Let us now investigate the contribution to the conduc- 
tivity arising from the regular part of the MT diagram. 
The operator for the electromagnetic response can be 
written as 



freg/ 



ilk 



__cos(g,d) 



2/i-£+r?(q) 



(3.35) 



As before, we can write this sum as a contour integral. 
Then, we have to perform the analytical continuation to 
real values of the external frequency ito^: iui^, lo. We 
expand the resulting expression in a Taylor-series up to 
the second order in uj. The static term is canceled by 
a similar term in (3(0). As a result we obtain for the 
operator of the electromagnetic response: 



r 



4ttT) 



coth 



n Qdn 



2T ML 



fj2 



(3.36) 



Performing the integration in Eq. (3.36) we obtain for 
the regular MT contribution to conductivity in the low 
temperature limit: 
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MT 
reg 
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i—1 



'7(q) 



(3.37) 



Let us estimate the integral in the right hand side of 
Eq. (3.37). At very low temperatures, T <C Agfy, and 
near the critical magnetic field He, h <C gS/Ao, it turns 
out that 



_5_ 

'A'o' 



(3.38) 



The regular MT part gives a negative contribution to the 
fluctuation conductivity and it is independent of temper- 
ature, so that even at zero temperature tr^fj remains 
finite and reduces the conductivity. We have shown be- 
fore that at very low temperature neither a^^ nor 
can change the monotonous increase of resistivity caused 
by the suppression of the density of states and the regular 
MT contribution. 



D. Diagrams No. 3, 4, 7, 8 and 9,10 

These diagrams arise when the DOS and MT dia- 
grams are averaged over the impurity positions. This 
averaging process results in the appearance of an ad- 
ditional Cooperon connecting two different grains with 
each other. Let us consider the diagram 9 in Fig. ^, the 
diagram 10 can be calculated in the same way. The ana- 
lytical expression of the operator of the electromagnetic 
response reads as follows 



i—l O/s 



(3.39) 



where K{ink,ci) is the propagator of the superconduct- 
ing fluctuations and B{iflk,ii-iJ^,q) corresponds to the 
contribution of the loop. A factor of 2 in Eq. (3.39) comes 
from the summation over spin indices and another factor 
of 2 originates from a similar diagram shown in Fig. |^. 
The analytical expression for the loop can be written as 



- ^ j -°('?'^)(^ j ^^"fa^) (2.)3 (2..or)3 ^g J 
xG{ien+u,Pi)G^{ien,P2)G{ink - ien,P2)C'^{ien,i^k - ien,q)C{ien+^,i^k ~ i£n,q), (3.40) 



where pi and P2 denote the momenta in the different 
granules and q',q" are quasimomenta. As before we can 
set rtk inside the Green's functions equal to zero. We 
can immediately see that after integration over q', q" the 
function B{i^lk, iuj^, q) is equal to zero and the contribu- 
tion to the fluctuation conductivity from this diagram 
vanishes. The same reason holds for the diagram 3 and 
4 in Fig. || which come from the averaging of the MT- 
diagram over impurity positions. Also the diagrams 7 
and 8 of Fig. |^ do not contribute to the fluctuation con- 
ductivity. For simplicity we present the diagram 7 once 
again in Fig. m 




P3 , 8+CO 

FIG. 4. DOS-type diagram with an additional impurity 
renormalization. 



One can easily see that the corresponding analytical 
expression of such a diagram contains a term of the fol- 
lowing form / ^|e|G'^(£,P2)G^(£,P2) = 0. Both Green 
functions have poles on the same side of the complex 
plane and therefore the integral is equal to zero. Thus, 
we may neglect diagrams 7 and 8. As we have shown, 
only the diagrams 1 (AL), 2 (MT), 5 and 6 (DOS) make 
a contribution to the fluctuation conductivity of granular 
metals. 



E. Final Formulae 

The calculations presented in the previous subsections 
show that in granular superconducting metals fluctua- 
tions make a considerable contribution to the conductiv- 
ity in the normal phase at low temperatures and strong 
magnetic field. The fluctuation conductivity of granu- 
lar metals is given by the DOS, Aslamazov-Larkin and 
Maki-Thompson contribution 

^fl ^ ^DOS ^ ^AL ^ ^MT _ (3 4^) 

Other contributions from the diagrams in Fig. ^ are equal 
to zero. At low temperatures, T <gC Aq??, the final result 
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for the total fluctuation conductivity of granular metals can be written as: 
I 



1 

2ir 



In 




47r ^_ 

9 \ 47rTaq>;^(q) 



sjqid) 



9ir 



^ ' cos(gi(i) In f ^ 



'7(q) 



(3.42) 



The Eq. (3.42) is the main result of our paper. For 
low temperatures, T <C Aqt) and strong magnetic field, 
H > He the main contribution comes from the fi rst and 
from the last terms in the brackets in Eq. ( 3.42| ). As a 
consequence, the superconducting fluctuation contribu- 
tion to the conductivity of granular metals is negative. 
The quantitative result of Eq. (3.42) is written for the 
region 1 <C (? <C In the next section we show that 
the correction to the conductivity due to weak localiza- 
tion effects can be neglected. 



Similar to the case of a homogeneous sample, the WL 
correction to conductivity has the negative sign. 

Now let us compare the WL correction with a^^^ . 
The ratio of these two quantities is given by 



5a 



WL 



^DOS 



-3/2 



Ao 

5 



3/2 



In" 



Ao 

95 



(4.3) 



For large g ^ 1, this ratio is small and therefore the WL 
correction is always smaller than the correction which 
arises from the suppression of the density of states. 



IV. WEAK LOCALIZATION CORRECTION TO 
CONDUCTIVITY OF GRANULAR METALS 



V. CONCLUSION 



The Weakly Localized Regime (WLR) is the regime 
where interference effects between different plane waves, 
being treated independently, start to play a role. The in- 
terference of plane waves leads to an increase of the prob- 
ability to find an electron at a certain place and this effect 
results in a reduction of the conductivity. The diageam 
corresponding to the WL correction to conductivitytj is 
shown in Fig. |^. 






= L 

FIG. 5. Weak localization correction to conductivity. The 
shaded block denotes the renormalized Cooperon. 

For the WL correction of granular metal^O we obtain 



^^geHd^Y. 



C(0,q) ^3 

cos(g,d)7— 3, (4.1) 



(27r)3 



here C(0, q) is the Cooperon taken at the frequency lo = 
and quasimomentum q. Using Eq. (3.4) and making the 
integration over quasimomentum q in Eq. (4.1) we obtain 
the final result for the WL correction to conductivity of 
granular metals 



5a 



WL 



-3/2 



Co 



Ao 



1/2 



(4.2) 



We have obtained the fluctuation conductivity of gran- 
ular metals at low temperatures and st rong magnetic 
field. Our main result is given by Eq. ( |3.42 ). To ob- 
tain this result we assumed that th e di mensionless con- 
ductance, 5, satisfies the inequality (|1 . l|) . Therefore, the 
granular structure of the metal is essential for our consid- 
eration. We have generalized the previous studieajQ to 
the entire region 1^5^ . This case is still different 
from the case of 2D homogeneous superconductorQ where 
the dimensionless conductance was assumed to be very 
large (of order kpl). 

One can see that a^^^ and cr^ej Si^^ ^ negative con- 
tribution to the fluctuation conductivity, whereas the AL 
and the anomalous MT contributions are positive. This 
leads to a competition between the positive and the neg- 
ative contributions. But as we have shown, at low tem- 
peratures, T <C Tc, and strong magnetic field, H ^ He, 
this negative contribution cannot be compensated by the 
positive contributions and therefore the entire fluctua- 
tion conductivity is negative. The situation holds even 
at r = where the AL and anomalous MT contribu- 
tions are equal to zero. 



Magnetic fieldh=(H-Hc2)/Hc2 



10 



FIG. 6. Schematic picture of the resistivity of a granu- 
lated superconductor at fixed temperature as a function of 
the magnetic field. 

Qualitatively results are depicted in Fig. ^, where the 
typical curve for the dependence of the fluctuation resis- 
tivity on the reduced magnetic field h at low temperature 
is presented. The curve reaches the value of the classi- 
cal resistivity Rq asymptotically onlyip extremely strong 
magnetic fields. It was shown in Ref.B that for the gran- 
ular metals the real transition into the superconducting 
state occurs not at but at a lower field H^^ , which is 
due to the electron motion over many grains. Thus, in or- 
der to take into account the macroscopic orbital electron 
motion we have to replace He by in our formulas. 
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APPENDIX A: 

In this appendix we carry out the analytical continu- 
ation in Eq. (3.11) and obtain the result for the opera- 
tor of electromagnetic response Q^{uj), Eq. (3.12). The 
lower limit in the sum over Qk in Eq- (3.11) depends 
on the external frequency lo^. Therefore the analytical 
continuation of this expression in the region of real fre- 
quencies is not correct by simply making the replacement 
iuji, u + lO"*", it can be achieved by a transformation 
of the sum into a contour integral 



T J2 Ditnk,iuj,,q)K{tnk,q)^ 

If z 
T~ / coth—D{z,iuj^,ci)K{z,q}dz, (Al) 



where the contours of integration Ci and C2 are shown in 
Fig. 1^. We should interpret the propagator of the super- 
conducting fluctuations as the retarded propagator 
when 5(2) > and as the advanced propagator when 
3(z) < 0. In the vicinity of the critical magnetic fleld, 
[H — Hc)/Hc ^ 1, we can expand the logarithm in the 



denominator of K and the retarded ptr advanced form of 
this quantity has the following forml 



here h = {H - Hc)/H, 



1 fr., , N 

- 2/iT-^+??(q) 



(A2) 




C2 

Im(z+iav )=0 



FIG. 7. Contours of integration Ci and C2, the arrows 
indicate the way going along the paths. 

The contribution to the integral from the large circle 
vanishes if we extend the contour to ioo, so that only 
the paths along the horizontal lines with Q(z) = and 
3(z + iLOjj) — remain. Thus for the right hand side in 
Eq. (Al) we obtain 



J coth— i?^(-^z,^c^,,q)if^(-^z,q)— (A3) 
f°° z 

\ coth— L»'^(-^z,^w^,,q)i^'^(-^z,q)- 
y-oo 2T 



—00 

+ 00 — iuj 



dz 

' A'Kl 



coth -^D^i-iz, iuj^, q)iv:^(-iz, q)^, 



where D^, are the retarded, advanced forms of the 
function D(iD,k, icj^, q). They are given by 



i?«-^(f7,c.,q) = - 



K^dT 



(A4) 



, , 1 2iLL! + iQ 



^ ' 2 ^ 47rr 



In the third integral in Eq. (A3) we make the substi- 
tution of variables z + iuj^ z and use the fact that 
coth (■^^^) = coth {^)- Now we can simply make the 
analytical continuation: iuj^ u. Finally we obtain for 
the sum over ilk in Eq. ( [Ai| ) the following result 



+00 



^ coth ^D'^Hn, q)i^«(-^^7, q) + /■ °° ^ Lth ^ - coth ^ ) Z?^(-^^], co, q)if ^(-z^l, q) 



(A5) 



Since we are interested in the dc conductivity it is sufficient to retain the linear term in u in Eq. (A5) and then we 
obtain: 
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+C50 +00 

i > luj — . / cotn— -w — + an -fi Q — / ■ 

4TT^dT3 J 2T^ \2 4ttT V ^ '^Uni 2tJ sinh' # 47ri 



(A6) 



The main contribution to the integral over in Eq. ( [A6| ) 
comes from small values of the frequency therefore we 
may put f2 = inside ipi^) the first integral and ex- 
tend the integrand by ifl/Ao + fj. One can easily see 
that this integral is logarithmically divergent and it has 
to be cut off at Qmax ^ Aq. In the second integral we 
make an expansion in Q, retaining only the first nonvan- 
ishing term and then we extend the resulting expression 
by — zfi/Ao + f]. As a result for t he op erator of electro- 
magnetic response we obtain Eq. ( 3.12| ). 
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